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ABSTRACT 


An  experimental  study  has  been  performed  of  the  unsteady  processes  in 
the  starting  period  of  a  supersonic  Ludwieg  tube,  a  device  which  operates 
like  an  intermittent  supersonic  wind  tunnel.  A  quick  opening  diaphragm  lo¬ 
cated  downstream  of  the  nozzle  initiates  the  flew.  Pressure  and  density 
measurements  are  made  in  a  variety  of  ways  in  Mach  number  1.67  and  3.0  noz¬ 
zles.  For  the  starting  conditions  treated,  supersonic  flow  is  established 
in  the  nozzle  without  producing  shock  waves.  Various  time  dependent  func¬ 
tions  are  observed  in  the  adjustment  of  gasdynamic  parameters  to  their  steady 
supersonic  values.  These  changes  of  pressure,  etc.,  include  undershoots, 
overshoots,  and  other  variations  of  the  final  steady-state  values.  Calcula¬ 
tions  based  an  an  assumed  zero- length  nozzle  do  not  adequately  predict  start¬ 
ing  times  and  pressures. 
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ILLUSTRATIONS 


Figure  1  Ludwieg  Tube:  Simplified  t  -  x  plot.  The  times 

schematically  indicated  are  t^,  the  total  starting 
time;  t-Q,  the  time  required  to  establish  sonic  con¬ 
ditions  at  the  nozzle's  throat;  tj,  the  total  time 
of  stable  supersonic  flow.  2 

Figure  2  Experimental  Apparatus.  (A)  High  Speed  Camera. 

(B)  Pressure  Transducer  (at  upstream  position) . 

(C)  Mach-Zehnder  Interferometer .  (D)  Side  View  of 

Test  Section.  Supply  tube:  5.30  in.  I.D.  Durrp  tube: 

3.76  in.  I.D.  Nozzle:  2"  x  2"  at  throat  (M  =  1.67  at 
exit) .  (E)  Collimating  Lens  and  Light  Source. 

(F)  Oscilloscope  to  record  pressure  traces.  5 

Figure  3  Sequential  Shadowgraphs  in  Starting  Stage.  Shock 

waves  are  produced  by  an  adhesive  strip  taped  to  the 

top  of  the  M  =  1.67  nozzle.  8 

Figure  4  Sample  Interferometer  Strip.  The  M  =  3.0  nozzle's 

throat  is  indicated  by  the  arrow.  Decreasing  densi¬ 
ties  are  indicated  by  the  motion  of  fringes  frcm  left 
to  right.  Flow  is  from  left  to  right.  9 

Figure  5  Sample  of  Density  Measurements.  Distances  are  measured 
in  inches  from  the  throat  of  the  nozzle.  Times  are 
measured  in  milliseconds,  with  t  =  0  corresponding  to 
first  appearance  of  the  expansion  fan  roughly  3.5  in. 
dcwnstream  of  the  throat.  10 

Figure  6  Static  Pressure  Measurements  in  the  M  =  1.67  Nozzle. 

The  pressure  gauge  is  being  self-triggered  in  these 
measurements;  the  signal  required  for  triggering  cor¬ 
responds  to  a  pressure  change  of  roughly  two  psi.  12 

Figure  7  Pressure  Measurements  in  the  Stagnation  Region: 

M  =  1.67.  (A)  P4  =  3.0,  P4/Px  =  7.7;  (B)  P4  =  3.0, 

P4/Pj  =  46;  (C)  P4  =  1.0,  P4/P!  =  7.7;  (D)  P4  =  3.7, 

P4/Pi  =  7.7.  The  time  scale  is  arbitrary;  the  pres¬ 
sure  trace  is  triggered  by  a  firing  pin  contact  sub¬ 
sequent  to  the  rupture  of  the  diaphragm.  P4  is  in 
atmospheres.  13 

Figure  8  Pressure  Measurements  in  the  Stagnation  Region: 

M  =  3.0.  (A)  P4  =  3.0,  P4/Pi  =  18;  (B)  P4  =  3.0, 

P4/Pj  =  46;  (C)  P4  =  2.0,  E^/Px  =  18;  (D)  P4  =  3.7, 

Rj/Px  =18.  The  pressure  gauge  triggering  here  is 

the  same  as  in  Figure  7.  P4  is  in  atmospheres.  14 
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Figure  9 


Starting  Times  and  Pressures  in  the  Stagnation 
Region:  Insensitivity  to  P4. 
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I 

INTRODUCTICN 


For  vehicles  which  fly  at  M  >  3  and  high  altitudes,  current  testing  faci¬ 
lities  can  simulate  Reynolds  numbers,  Mach  numbers  and  flight  environments 
with  adequate  success  so  that  flight  performance  can  be  predicted  fran  model 
and  full  scale  experiments  (e.g.  Ref.  1) .  However,  in  the  Mach  number  range 
0.8  <  M  <  3,  this  is  less  true.  The  high  Reynolds  numbers  experienced  by  air¬ 
craft  and  space  vehicles  cannot  be  satisfactorily  reproduced  in  the  laboratory* 
For  example,  at  M  ~  3.5  the  Saturn  V  flies  at  Rg  -  2  x  106  whereas  at  M  ■  1.5 
the  Reynolds  number  is  as  high  as  10 7  (based  on  vehicle  length) .  Throughout 
this  Mach  number  range,  testing  facilities  currently  available  produce  Reynolds 
numbers  of  only  roughly  2  x  10  .  Since  Reynolds  number  simulation  under  these 
circumstances  in  a  crucial  aspect  of  vehicle  design,  considerable  interest  has 
developed  in  attempts  to  produce  testing  facilities  with  high  stagnation  pres¬ 
sure,  long  test  time  at  transonic  and  moderately  supersonic  speeds. 

Several  years  ago,  H.  Ludwieg  suggested  in  a  different  context  a  device 
which  now  appears  to  be  remarkably  well  suited  for  this  kind  of  application 
(2) .  The  device  consisted  of  a  long  cylindrical  tube  as  a  container  for  the 
compressed  air.  One  end  of  the  tube  is  closed;  the  other  end  contained  a 
supersonic  nozzle,  a  test  section  and  a  quick  opening  valve  which  is  opened 
to  the  atmosphere.  In  the  Yale  University  Ludwieg  tube  (3)  a  conventional 
shock  tube  configuration  is  modified  by  the  insertion  of  a  supersonic  nozzle 
into  the  section  upstream  of  the  diaphragm  (the  high  pressure  section) . 

After  the  diaphragm  breaks  a  shock  wave  and  a  contact  surface  travel  down¬ 
stream  and  the  head  of  the  expansion  wave  moves  upstream  through  the  super¬ 
sonic  nozzle.  After  the  remaining  part  of  the  expansion  wave  is  swept  back 
downstream,  stable  conditions  of  supersonic  flew  are  maintained  in  the  nozzle 
until  the  reflected  expansion  wave  returns  to  the  throat.  The  operating  stages 
for  this  intermittent  tube-wind- tunnel  (also  knewn  as  the  "Ludwieg  tube")  cure 
indicated  in  Figure  1,  a  t  -  x  diagram  with  a  sketch  of  the  tube.  When  the 
diaphragm  is  ruptured,  a  centered  expansion  wave  is  formed  and  moves  upstream 
into  the  high  pressure  tube  where  the  gas  is  at  the  initial  conditions  denoted 
by  4.  Part  of  it  passes  through  the  nozzle  and  travels  to  the  end  of  the  high 
pressure  tube  where  it  is  reflected  and  returns  to  the  nozzle.  This  marks  the 
end  of  the  first  steady  flow  period.  When  the  diaphragm  breaks,  a  shock  wave 
is  formed  and  travels  downstream  followed  by  the  contact  surface  into  the  dump 
tube  where  the  initial  conditions  of  the  gas  are  denoted  by  1.  The  times 
shown  are  tj_,  the  total  starting  time;  tj_]_,  the  time  required  to  establish 
sonic  conditions  at  the  nozzle's  throat;  and  bp,  the  total  time  of  stable 
supersonic  flow.  In  the  literature,  one  can  find  derivations  of  the  formulas 
relating  the  nozzle  parameters  to  the  initial  conditions  in  the  tube.  In  these 
analyses,  it  is  generally  assumed  that  the  nozzle  has  zero  length,  that  the 
flow  is  one-dimensional  and  inviscid  and  that  one  is  dealing  with  a  perfect 
gas  as  shown  e.g.  by  Cable  &  Cox  (4)  and  Becker  (5) .  With  this  device,- high 
stagnation  pressures  can  be  maintained  in  the  high  pressure  side  and  steady 
supersonic  flow  of  relatively  short  but  useful  duration  can  be  established  in 
the  test  section.  Furthermore,  a  practical  advantage  with  respect  to  high 
pressure  blow-down  tunnels  arises  fran  the  fact  that  no  valves  are  required 
to  maintain  steady  nozzle  supply  conditions .  For  these  reasons,  interest  new 
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Figure  1  Ludwieg  Tube:  Simplified  t  -  x  plot.  The  times  schematically  indicated 

are  tj_,  the  total  starting  time;  tn,  the  time  required  to  establish  sonic 
conditions  at  the  nozzle's  throat;  bj,  the  total  time  of  stable  supersonic 
f: 
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exists  in  the  Ludwieg  tube  and  in  the  nature  of  its  possible  application. 

With  regard  to  the  starting  stage,  that  is,  the  elapsed  time  between  the 
rupture  of  the  diaphragm  or  the  opening  of  a  valve  and  the  establishment  of 
steady  supersonic  flow  in  the  nozzle,  sene  uncertainties  have  persisted. 
First,  there  is  the  question  of  whether  or  not  starting  shock  waves  are  pro¬ 
duced  after  sonic  conditions  are  obtained  at  the  throat  of  the  nozzle.  Such 
shocks  may  originate  at  the  throat  and  travel  downstream  through  the  nozzle 
prior  to  the  steady  flow  stage  and  may  have  adverse  effects  on  the  model.  In 
addition,  there  is  the  question  of  how,  in  fact,  the  steady  state  parameters 
in  the  nozzle  are  acquired. 

In  this  paper  we  shall  present  further  initial  results  from  our  obser¬ 
vations  of  starting  processes  in  a  supersonic  Ludwieg  tube*.  Our  approach 
has  been  phenomenological;  we  are  concerned  at  this  point  only  with  experi¬ 
mentally  determining  what  seen  to  us  to  be  the  gross  qualitative  features  of 
this  starting  phase.  We  are  interested  in  observing  the  sensitivity  of  the 
starting  processes  to  changes  in  nozzle  pressure,  tube  pressure  ratio  across 
the  unruptured  diaphragm,  and  nozzle  Mach  number. 


*Preliminary  results  have  been  given  in  Reference  6. 
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II 

EXPERIMENTAL  APPARATUS 

Hie  Yale  University  Ludwieg  tube  (3)  is  sketched  in  Figure  2.  The  over¬ 
all  length  is  26  ft,  and  a  diaphragm  can  be  placed  on  either  side  of  the  test 
section*.  The  two  dimensional  wedge-nozzle  used  in  most  of  these  studies  has 
an  exit  Mach  number  of  1.67  and  a  throat  to  exit  distance  of  six  inches  with 
an  expansion  half -angle  of  three  degrees.  It  can  be  placed  at  various  posi¬ 
tions  in  the  1.5  ft  test  section  so  as  to  vary  the  field  of  view.  Two  sec¬ 
tions  immediately  upstream  and  downstream  of  the  test  section  respectively 
produce  a  smooth  transition  frcm  the  circular  supply  and  dump  tubes  bo  the 
rectangular  test  section.  In  all  experiments,  the  operating  gas  is  dry  nitro¬ 
gen  or  dry  air. 

Three  diagnostic  techniques  have  been  used  in  our  investigations: 

(1)  Timed  shadowgraphs  are  made  during  the  starting  stage.  High  speed 
movie  shadowgraphs  with  a  0.25  msec  frame  to  frame  sampling  rate  are  taken 
losing  a  Fas  tax  camera.  Spark  shadowgraphs  are  taken  at  0.1  msec  intervals 
using  a  1  ysec  10,000  volt  spark  source. 

(2)  Slit  streak-interferometry  is  used  to  obtain  density  measurements 
as  a  function  of  time  on  the  nozzle  center-line.  A  horizontal  slit  is  placed 
between  the  light  source  and  the  collimating  lens  so  that  its  image  appears  at 
the  test  section  with  a  width  of  2  mm.  This,  together  with  the  equivalent 
film  speed  of  300  m/sec,  gives  a  time  resolution  of  roughly  7  ysec.  The  Zeiss 
Mach-Zehnder  interferometer**  (plate  size;  4  1/8"  x  7")  is  adjusted  so  that 
the  image  of  vertical  interference  fringes  appears  at  the  center  of  the  test 
section.  With  the  rotating  prism  in  the  Fas  tax  movie  camera  removed  and  with 
tiie  film  moving  in  a  vertical  plane,  continuous  interferograms  are  obtained 
which  allcw  one  to  follow  the  movement  of  each  fringe,  thus  permitting  a 
determination  of  the  change  in  density  with  time  at  any  position  along  the 
centerline  of  the  nozzle.  Shock  waves  would  be  noticeable  as  discontinuous 
jumps  in  the  otherwise  smooth  increase  and  decrease  in  density  that  can  be 
seen  according  to  the  direction  of  motion  of  the  fringes.  The  system  is  found 
to  be  sensitive  to  density  changes  of  less  than  1.0%. 

(3)  Static  pressure  measurements  have  been  made  in  the  supply  tube  and 
at  various  locations  in  the  nozzle  using  a  calibrated  high  speed  3  ysec  rise- 
time  quartz  pressure-transducer  (Kistler,  Model  #606L) . 

For  a  more  detailed  description  of  the  tube  wind  tunnel  and  auxiliary 
equipment,  see  Appendix  I. 


‘Downstream  diaphragm  used  in  all  experiments. 

**We  are  grateful  to  the  U.S.  Naval  Ordnance  Laboratory  for  loaning  us  this 
instrument. 
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Figure  2  Experimental  Apparatus..  (A)  High  Speed  Camera.  (B)  Pressure  Transducer'  (at  upstream 
position) .  (C)  Mach-Zehnder  Interferometer.  ■  (D)  Side  View  of  Test  Section.  Supply 
Tube:  5.30  in.  I.D.  Dump  tube:  3.76  in.  I.  D.  Nozzle:  2"  x  2"  at  throat  ;(M  =  1.67 
.  at- exit).  ,  (E)  Collimating  Lens  ■  and  Light' Source'.  (F)  Oscilloscope  to  record  pressure 
-  .  traces.  :  .  '  -  ~ 
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III 

EXPERIMENTAL  RESULTS 


Spark  shadowgraph  pictures  are  taken  during  steady  flew  with  a  model 
inserted  in  the  test  section  producing  bow  shock  waves  to  confirm  that  super¬ 
sonic  flow  is  established  (see  Appendix  II) .  Movie  shadowgraphs  are  taken 
with  fiducial  marks  on  the  test  section  window  and  with  adhesive  strips  along 
the  diverging  section  of  the  nozzle  and  perpendicular  to  the  flow.  The  rela¬ 
tive  position  of  the  disturbances  frem  the  adhesive  strips  with  the  marks  on 
the  window  indicates  when  supersonic  flew  begins  and  reaches  steady  state. 
Finally,  movie  shadowgraphs  and  spark  shadowgraphs  of  the  flew  in  the  nozzle 
permit  a  search  to  be  made  for  starting  shocks  in  the  starting  stage.  Table 
I  lists  the  operating  pressures  for  the  experiments  with  movie  shadowgraphs. 
When  mare  than  one  movie  is  indicated,  this  means  that  different  regions  in 
the  nozzle  have  been  sampled  or  that  models  have  been  inserted  into  the  test 
section.  In  addition,  a  set  of  spark  shadewgraphs  has  been  taken  during  the 
starting  stage  for  P4  =  1.0  atm  and  P4/P1  =  7.7  with  a  sampling  rate  of  better 
than  one  every  0.1  msec. 

In  none  of  these  experiments  have  starting  shock  waves  been  observed. 

This  is  illustrated  in  Figure  3  vdiere  P4  =  3.0  atm  and  P4/P1  =  23  for  the 
M  =  1.67  nozzle.  Here  one  sees  supersonic  flew  being  established  in  the  noz¬ 
zle  with  no  waves  other  than  the  disturbances  from  the  adhesive  strips  on  the 
nozzle  block.  The  supersonic  flew  appears  to  develop  smoothly  at  all  points 
in  the  nozzle.  Similar  results  are  obtained  for  all  cases  studied  by  us.  In 
Figure  4,  one  sees  an  example  of  the  movie  interferogram  by  which  one  can  look 
for  starting  shocks  with  even  greater  reliability.  Furthermore,  density  his¬ 
tories  at  fixed  nozzle  position  and  density  profiles  at  specific  times  in  the 
starting  stage  can  be  determined.  The  experiments  during  which  movie  inter- 
ferograms  were  made  are  also  indicated  in  Table  I.  The  range  of  upstream 
pressures  and  of  pressure  ratios  is  nearly  the  same  as  that  for  the  movie 
shadewgraphs.  In  Table  I  one  sees  that  neither  the  interferograms  nor  the 
shadewgraphs  shew  starting  shocks.  With  a  total  of  seventeen  experiments  per¬ 
formed,  including  five  experiments  at  P4  =  3.0  atm  and  P4/P1  =46,  one  finds 
no  indication  in  these  data  of  a  starting  shock  prior  to  the  establishment  of 
supersonic  flew. 

Examples  of  density  measurements  are  given  in  Figure  5.  The  throat  of 
the  nozzle  is  at  X  =  0  and  the  time  scale  is  arbitrarily  chosen  so  that 
t  =  0  corresponds  to  the  first  appearance  of  the  head  of  the  expansion  fan 
roughly  3.5  in.  dewnstream  of  the  throat.  One  can  see  the  density  drop  to 
values  that  are  at  first  below  those  of  its  final  steady  state.  This 
"undershoot"  is  about  4%  at  X  =  0;  however,  its  magnitude  decreases  with  in¬ 
creasing  dewnstream  distance.  One  also  notices  that  the  undershoot  occurs 
in  about  2  msec  at  the  throat  and  that  it  precedes  the  appearance  of  an  under¬ 
shoot  at  the  dewnstream  distance  for  X  >  0.  Density  profiles  are  also  shown 
in  Figure  5  for  this  case  of  P4  =  3.0  atm  and  Pj  =  0.40  atm.  At  all  times, 
the  densities  show  a  decrease  with  increasing  dewnstream  distance.  Here,  as 
before,  there  is  no  indication  of  starting  shock  waves.  Comparisons  of  the 
density  measurements  in  Figure  5  with  density  measurements  at  other  values  of 
P4  (with  P4/P1  =  7.7  =  const.)  have  given  the  following  results:  The  values 
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Table  I.  Summary  of  Experiments  M  =  1.67 


(atm) 

P4A1 

No.  of  Movie- 
Shadcwgraphs 

No.  of  Movie- 
Interferograms 

Starting  Shocks 
Formed? 

3.0 

46 

3 

2 

No 

3.0 

23 

1 

— 

No 

3.0 

12 

1 

— 

No 

3.0 

7.7 

2 

2 

No 
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t  =  t  +  1.4  msec  t  =  t  +1.6  msec  t 
o  o 


t  +1.8  msec 
o 


Figure  3  Sequential  Shadowgraphs  in  Starting  Stage. 

Shock  waves  are  produced  by  an  adhesive  strip 
taped  to  the  top  of  the  M  =  1.67  nozzle. 
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Figure  5  Sample  of  Density  Measurements.  Distances  are  measured  in  inches  from  the  throat 
of  the  nozzle.  Times  are  measured  in  milliseconds,  with  t  =  0  corresponding  to 
first  appearance  of  the  expansion  fan  roughly  3.5  in.  downstream  of  the  throat. 
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of  p/o 4  obtained  at  the  throat  are  independent  of  P4;  hence  the  percentage 
undershoot  is  also  independent  of  P4.  The  undershoot  times  are  independent  of 
P4.  Further,  the  density  profiles  always  show  a  monotonic  decrease  with  in¬ 
creasing  downstream  distance.  Generally,  we  have  found  that  the  data  in 
Figure  5  are  representative  of  all  density  measurements  performed. 

Direct  measurements  of  static  pressure  histories  in  the  nozzle  provide 
additional  checks  on  the  existence  or  absence  of  starting  shock  waves.  In 
addition,  comparisons  can  be  made  with  the  density  histories  mentioned  above. 
Pressure  as  a  function  of  time  at  three  different  locations  in  the  M  =  1.67 
nozzle  is  given  in  Figure  6.  The  starting  point  of  the  time  scale  is  arbi¬ 
trary.  At  X  =  0  in.  the  pressure  drops  initially  below  its  final  value  while 
at  X  =  5.2  in.  this  undershoot  is  considerably  diminished;  outside  the  noz¬ 
zle  at  X  =  6.3  in.  this  behavior  has  practically  disappeared.  These  results 
are  consistent  with  the  density  measurements  shewn  in  Figure  5  where  the  re¬ 
lative  undershoot  also  decreases  with  increasing  downstream  distance.  Density 
measurements  made  at  the  nozzle's  exit  show  the  same  qualitative  trends  as  the 
pressure  measurements  in  Figure  6.  Finally,  Figure  6  provides  additional 
evidence  for  the  smoothness  and  shock-free  nature  of  the  establishment  of 
supersonic  flow  in  the  nozzle. 

Measurements  of  static  pressure  histories  in  the  stagnation  region  have 
been  obtained  for  all  the  experiments  given  in  Table  I.  Samples  of  these  re¬ 
sults  are  given  in  Figure  7  for  the  M  =  1.67  nozzle*.  An  undershoot  of  the 
steady  state  values  is  also  seen  as  a  persistent  feature.  In  addition,  one 
finds  that  the  undershoot  seems  to  be  followed  by  an  "overshoot"  as  the  ad¬ 
justment  to  steady  flew  values  is  taking  place.  By  acmparing  the  two  results 
for  P4  =  3.0  atm  and  different  Pi  values,  one  can  see  that  Px,  the  initial 
downstream  tube  pressure,  is  not  the  determining  factor  as  far  as  the  pressure 
history  is  concerned,  once  Pi  is  low  enough  to  provide  a  sufficient  pressure 
ratio  for  supersonic  flow  to  be  established.  Changes  in  P4  do  cause  changes 
in  the  magnitude  of  the  pressure  shifts  during  starting.  However,  the  times 
required  for  the  pressure  to  drop  below  its  final  values  seem  to  be  indepen¬ 
dent  of  both  P4  and  P4/Pi  within  the  limits  stated. 

For  purposes  of  comparison,  static  pressure  measurements  have  also  been 
made  in  the  stagnation  region  of  a  taro  dimensional  Mach  number  3.0  nozzle. 
Samples  of  the  data  obtained  are  given  in  Figure  8.  Qualitatively,  the  re¬ 
sults  are  similar  to  those  observed  at  the  lower  Mach  number. 

For  a  more  detailed  description  of  the -experiments  and  further  examples 
of  the  results,  see  Appendix  II. 


*Static  and  stagnation  pressure  measurements  in  the  stagnation  region  are  also 
reported  in  References  7  and  3.  However,  these  experiments  were  designed  to 
determine  steady  state  values  and  gave  little  information  on  the  starting 
processes. 
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(2msec/div) 


Figure  6  Static  Pressure  Measuranents  in  the  M  =  1.67  Nozzle.  The 
pressure  gauge  is  being  self -triggered  in  these  measure¬ 
ments;  the  signal  required  for  triggering  corresponds  to  a 
pressure  change  of  roughly  two  psi. 
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M  =  l.67 
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Fiaure  7  Pressure  Measurements  in  the  Stagnation  Region:  M  -  1*67.  (A)  _p4  •  • 

p./p.  =  7.7;  (B)  P4  =  3.0,  P4/P1  =  46;  (C)  P4  =  1.0,  P4/Pi  -  7.7; 

(D)  P4  =  3.7,  P4/Pi  =7.7.  The  time  scale  is  arbitrary;  the  pressure 
trace  is  triggered  by  a  firing  pin  contact  subsequent  to  the  rupture  of 
the  diaphragm.  P4  is  in  atmospheres. 
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Figure  8  PressuLe  Measurements  in  me  Stagnation  Region:  M  =  3.0.  (A)  P4  =  3.0, 

P4/P!  =  18;  (B)  P4  =  3.0,  P4/Pi  =  46;  (C)  P4  =  2.0,  P4/Pi  =  18; 

(D)  P4  =  ‘3.7,  P4/P;  =  18.  The  pressure  gauge  triggering  here  is  the  same 
as  in  Figure  7.  P4  is  in  atmospheres. 
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IV 

DISCUSSION  OF  RESULTS 


For  the  flew  conditions  of  this  work,  it  has  been  shown  that  starting 
shock  waves  analogous  to  those  which  ordinarily  appear  in  a  conventional 
blowdown  wind  tunnel  are  not  produced.  The  absence  of  starting  shocks  is 
important  because  it  reduces  the  transient  forces  on  a  model  during  tunnel 
starting  and  thereby  decreases  the  chance  of  model  breakage  and  lowers  the 
design  requirements  of  the  model  and  its  support.  The  M  =  1.67  nozzle  used 
in  this  work  is  a  wedge  type  nozzle  with  a  6  in.  radius  arc  at  the  throat  and 
a  straight  3°  half-angle  diverging  section  and  the  M  =  3.0  nozzle  is  a  contin¬ 
uously  expanding  nozzle  with  a  cubic  profile,  as  oempared  to  nozzles  designed 
by  using  the  method  of  characteristics  which  produce  parallel  supersonic  flow 
at  the  exit.  This  nozzle  design  may  effect  the  results.  We  recall  that  the 
diaphragms  used  in  these  experiments  require  less  than  0.1  msec  to  burst.  By 
contrast  roughly  1  msec  is  required  for  the  head  of  the  expansion  wave  to  tra¬ 
verse  the  nozzle.  This  time  is  called  the  characteristic  nozzle  time.  Thus, 
the  diaphragm  rupturing  process  is  fast  relative  to  the  characteristic  nozzle 
time.  On  the  other  hand,  the  valve  opening  process  in  a  blowdown  wind  tunnel 
is  typically  slew  relative  to  characteristic  nozzle  times.  It  is  therefore 
possible  that  a  Ludwieg  tube  lasing  slew  opening  diaphragms  or  valves  and  a 
short  nozzle  would  produce  starting  shock  waves.  For  the  M  =  1.67  nozzle, 
boundary  layer  effects  are  negligible  in  the  experiments  using  supply  pres¬ 
sures  between  1.0  and  3.7  atmospheres (5) .  This  conclusion  follows  from  the 
density  measurements  obtained  during  the  steady  supersonic  flow  in  the  nozzle. 
It  was  found  that  the  Mach  number  determined  from  the  density  ratio  was 
practically  equal  to  that  predicted  from  the  geometrical  area  ratio. 

However,  the  theory  for  the  Ludwieg  tube's  operating  parameters  for  zero 
nozzle  length  (e.g.  Reference  4)  is  not  entirely  adequate.  For  example,  it 
predicts  for  M  =  1.67  that  the  steady  state  static  pressure  should  be  86%  of 
the  supply  pressure  (Pst/P4  =  0*86)  while  the  measured  result  is  Pst/P4  “ 

0.83  ±  .01.  Hie  predicted  duration  of  the  stead/  flow  at  the  stagnation 
regions  pressure  station  is  17  msec  The  corresponding  measured  time  between 
the  appearance  of  the  head  of  expansion  wave  and  its  reflection  which  ter¬ 
minates  the  steady  flow  agrees  with  the  prediction.  However,  due  to  the 
starting  process  whose  duration  varies  from  5  to  7  msec  (depending  on  Pi,}, 
the  actual  time  of  stead/  flow  is  only  about  12  msec. 

It  was  found  that  we  could  make  some  correct  predictions  by  using  a 
simplified  one-dimensional  treatment  (see  Appendix  III) .  In  these  calcula¬ 
tions,  we  use  the  method  of  characteristics  applied  to  our  tube  configuration 
for  unsteady  one-dimensional  flow  in  the  M  =  1.67  nozzle.  It  was  assumed 
that  the  flew  was  initiated  by  an  expansion  fan  at  the  diaphragm  position. 

The  gas  was  again  assumed  to  be  perfect  and  inviscid.  It  was  then  found 
theoretically  that  1.6  msec  after  the  arrival  of  the  head  of  the  expansion 
wave  at  the  throat  of  the  nozzle,  the  density  at  the  throat  (X  =  0)  is 
p/p 4  =  0.51.  This  can  be  compared  with  the  measured  undershoot  time  and 
density  of  t  =  1.9  msec  and  p/p 4  =0.52  respectively.  Both  theory  and  experi¬ 
ment  give  the  steady  state  values  p/p 4  =  0.57  and  therefore  the  initially 
lower  density  value  is  verified  theoretically. 
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In  Figure  9,  finally,  two  quantitative  aspects  of  the  undershoot  in  the 
static  pressure  measurements  in  the  stagnation  region  are  shewn.  First  of 
all,  it  is  suggested  in  our  data  that  the  ratio  of  minimum  pressure  to  the 
steady  state  pressure  is  independent  of  the  upstream  starting  pressure.  For 
a  given  nozzle,  we  shew  in  Figure  9  data  at  a  single  fixed  pressure  ratio. 
However,  the  measurements  at  pressure  ratios  of  46,  23,  and  12  for  M  =  1.67 
and  of  46  for  M  =  3.0  show  these  results  to  be  indeed  independent  of  pressure 
ratio.  Secondly,  the  experiments  suggest  that  the  undershoot  time  and  the 
overshoot  time  are  independent  of  the  starting  pressure.  This  is  also  shown 
in  Figure  9.  The  time  between  the  onset  and  the  termination  of  the  expansion 
phase  in  the  stagnation  region  is  observed  to  be  insensitive  to  changes  in 
the  upstream  supply  pressure.  Nonetheless,  the  total  starting  time,  not 
shown  in  these  presentations,  does  seem  to  be  sensitive  to  Pi*. 
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Figure  9  Starting  Times  and  Pressures  in  the  Stagnation  Region: 
Insensitivity  to  Pu. 
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V 

SUM^RY 


From  our  experimental  and  theoretical  studies  of  the  starting  processes 
for  the  Ludwieg  tube  wind-tunnel  with  two  nozzles  of  M  =  1.67  and  M  =  3.0 
respectively,  for  the  range  of  conditions  treated  in  our  experiments,  we  find: 
(1)  Steady  supersonic  flew  can  be  established  smoothly  throughout  our  type 
of  continuously  diverging  nozzle  past  the  throat,  without  the  formation  of 
starting  shock  waves.  When  sonic  conditions  have  moved  to  the  throat  of  the 
nozzle,  supersonic  flow  has  already  begun  in  the  downstream  portion  of  the 
nozzle.  These  results  appear  to  be  independent  of  the  supply  pressure  P4  and 
of  the  starting  pressure  ratio  P4/P1,  once  this  ratio  is  sufficiently  large 
to  produce  supersonic  flew.  (2)  Static  pressures  and  densities  are  found  to 
first  undershoot  then  adjust  to  their  steady  supersonic  flow  values.  The 
values  of  the  pressure  and  density  drop  below  the  final  values  decrease  with 
increasing  distance  from  the  nozzle's  throat.  For  a  fixed  Mach  number  the 
ratios  of  minimum  pressure  and  density  to  supply  tube  pressure  and  density 
respectively  and  the  ratios  of  time  to  minimum  pressure  to  the  total  starting 
time  are  independent  of  P4  and  P4/P1.  (3)  For  fixed  nozzle  Mach  numbers, 

the  ratios  of  steady  state  pressure  to  supply  tube  pressure  are  independent 
of  P4  and  P4/P1.  However,  the  total  starting  time  is  found  to  be  dependent 
on  P4  but  not  on  P4/P1.  (4)  A  change  in  nozzle  Mach  number  changes  under¬ 

shoot  and  steady  state  pressures,  densities  and  time. 

These  results  suggest  that  the  deleterious  effects  arising  from  the  in¬ 
teraction  of  starting  shock  waves  with  test  models  can  be  avoided  in  Ludwieg 
tube  wind-tunnels  where  use  is  made  of  fast  opening  diaphragms  or  valves. 

The  observed  variations  in  pressures  and  densities  indicate  important  in¬ 
adequacies  in  a  simplified  one-dimensional  treatment  of  the  gas  dynamic  pro¬ 
cesses  and  show  that  more  sophisticated  theoretical  techniques  are  needed  if 
correct  predictions  concerning  the  starting  stage  phenomena  are  to  be  obtained. 
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DETAILED  DESCRIPTION  OF  LUDWIEG  TUBE  WIND  TUNNEL 
AND  AUXILIARY  EQUIPMENT 


A  detailed  side  view  of  the  tunnel  is  shewn  in  Figure  1-1.  The  overall 
length  of  the  intermittent  tunnel  is  26  ft.  and  it  is  fabricated  entirely  out 
of  stainless  steel  304.  The  high  pressure  end  of  the  tunnel  is  made  from  a 
10  ft.  length  of  seamless  tube  with  5.295  in.  I.D.  and  0.134  in.  wall  thick¬ 
ness.  A  transition  piece  made  out  of  1/8  in.  plate  and  1  ft.  in  length  con¬ 
nects  the  tube  to  a  rectangular  nozzle  section  1-1/2  ft.  in  length  and  con¬ 
tains  a  diaphragm  rupturing  pin.  The  high  pressure  end  of  the  tunnel  and  the 
transition  piece  are  mounted  on  linear  ball  bushings  which  make  it  possible 
by  a  simple  hand  push  to  move  the  entire  tube  and  transition  piece  4  in. 
away  from  the  nozzle  section  along  the  tunnel  centerline.  This  permits  in¬ 
sertion  and  removal  of  diaphragms  for  experiments  with  upstream  diaphragm 
location.  Screens  may  be  introduced  at  this  point.  The  nozzle  section  fol¬ 
lows  and  will  be  described  in  greater  detail  be  lew.  For  experiments  with 
diaphragm  downstream  of  nozzle  and  test  section  the  diaphragm  is  placed  be¬ 
tween  the  nozzle  section  and  another  1  ft.  transition  section  which  contains 
the  diaphragm  rupturing  pin  and  connects  the  nozzle  section  to  the  low  pres¬ 
sure  or  dump  tube.  The  dump  tube  is  made  of  seamless  tube  3.760  in.  I.D.  and 
0.120  in.  wall  thickness  and  consists  of  two  sections.  The  first  section  is 
2-1/2  ft.  long  and  is  connected  to  the  transition  piece.  The  other  section 
is  10  ft.  in  length.  A  2  in.  thick  insert,  the  size  of  a  4  in.  pipe  flange, 
is  located  between  these  two  pipes.  The  insert  contains  two  quartz  windows, 
3/4  in.  diameter  and  1  in.  thick,  and  permits  a  beam  of  light  to  pass  through 
the  tube.  This  arrangement  acts  as  a  shock  detection  station  and  can  be  used 
to  trigger  the  electronic  instrumentation.  The  individual  sections  of  the 
tube  are  sealed  together  by  inserts  containing  0-rings  placed  between  the 
flanges . 

The  downstream  transition  section  and  dump  tube  assembly  are  mounted  on 
linear  ball  bushings  which  make  it  possible  to  move  the  entire  assemble  ap¬ 
proximately  1  ft.  away  from  the  nozzle  section  along  the  tunnel  centerline. 
This  permits  easy  access  to  the  test  section  as  well  as  insertion  and  removal 
of  diaphragms.  The  nozzle  section,  diaphragm  and  downstream  transition  sec¬ 
tion  are  held  together  by  a  wedge-shape  clamping  device.  The  clamping  device 
permits  quick  connection  and  disconnection  of  the  two  parts  of  the  tunnel. 

The  entire  tunnel  assembly  is  mounted  on  a  4  in.  I-beam  and  supported  by 
several  4  in.  diameter  pipes.  This  arrangement  is  rigid  and  permits  easy 
access  to  any  part  of  the  tube.  A  3  ft.  piece  of  the  supporting  I-beam  for 
the  test  section  is  removable.  The  nozzle  section  and  this  portion  of  the 
I-beam  may  be  lifted  and  the  Mach-Zehnder  interferometer  (plate  size;  4  1/8"  x 
7") ,  manufactured  by  Zeiss  Company,  positioned  in  place. 

The  nozzle  section  is  rectangular  with  2  in.  by  5  in.  inside  dimensions 
and  1-1/2  ft.  in  length.  This  section  contains  the  nozzle  blocks,  test  sec¬ 
tion  and  model  support.  The  side  walls  extend  to  practically  the  entire 
length  of  the  section  and  are  sealed  to  it  by  means  of  0-rings.  The  nozzle 
blocks  are  2  in.  in  width  and  rest  on  the  top  and  bottom  of  the  rectangular 
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Pressure  tube:  5.295  in.  I.D.,  0.134  in.  wall,  stainless  steel  304  seamless  tube. 
Dump  tube.  3.760  in.  I.D.,  0.120  in.  wall,  stainless  steel  304  seamless  tube. 
Nozzle:  two  dimensional,  2x2  in.  throat. 


Figure  1.1  Sketch  of  the  intermittent  tube  wind  tunnel 
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nozzle  section.  They  are  held  in  place  by  screws  and  are  sealed  by  0-rings 
concentric  with  the  screws.  The  side  walls  may  be  easily  removed  to  exchange 
the  nozzle  blocks  which  have  been  machined  out  of  Plexiglass .  The  windows 
are  mounted  and  sealed  in  metal  frames  which  in  turn  are  sealed  by  means  of 
0-rings  into  the  side  walls.  Windows  are  held  flush  with  the  inside  of  the 
nozzle  walls  by  threaded  retaining  rings.  The  windows  are  4  in.  diameter, 

1  in.  thick,  and  made  of  Schott  BK-7  glass  ground  to  1/S  wave  length  with  a 
wedge  angle  no  greater  than  20  min.  The  location  of  the  windows  is  5  in.  up¬ 
stream  from  the  diaphragm.  The  model  support  is  mounted  in  the  nozzle  sec¬ 
tion  near  the  diaphragm.  Models  extend  into  the  nozzle  section,  as  viewed 
through  the  window,  and  are  easily  changed  by  reaching  into  the  test  section 
while  the  diaphragm  is  being  changed.  For  a  more  detailed  description  of  the 
equipment  see  Reference  3. 

Instrumentation  of  the  Tunnel 

The  tunnel  is  instrumented  with  a  pressure  measuring  system,  a  photo¬ 
graphic  system,  a  Mach-Zehnder  Interferometer,  and  a  supporting  electronic 
system.  See  Figure  2. 

Transient  pressure  recordings  are  made  with  a  piezo-electric  quartz- 
crystal  pressure  transducer  (Kistler  Model  606L) .  The  transducer  may  be 
mounted  in  an  insert  between  the  pressure  tube  and  the  upstream  transition 
section,  at  the  pressure  tube  end  plate,  and  also  in  the  nozzle  side  wall. 

The  side  wall  location  is  5.5  in.  upstream  from  the  diaphragm  and  the  trans¬ 
ducer  is  mounted  in  a  blank  which  replaces  the  window  when  pressure  readings 
are  taken.  The  output  of  the  transducer  is  amplified  by  a  Kistler  Universal 
Dial-Gain  Charge  Amplifier  Model  504  and  it  is  recorded  on  an  oscilloscope. 

The  rise  time  is  given  to  be  3  microseconds.  Simultaneous  pressure  record¬ 
ings  at  two  different  locations  are  obtained  by  using  two  pressure  transducers 
and  amplifiers.  Pressure  recordings  are  made  along  the  nozzle  centerline  at 
five  positions  each  3/4  in.  apart  by  using  a  special  mounting  plate  in  the 
nozzle  side  wall. 

The  Mach-Zehnder  Interferometer  is  located  at  the  test  section  and  is 
used  for  density  measurements .  This  is  done  by  streak  and  still  interfero- 
grams.  It  may  also  be  modified  so  that  shadowgraph  pictures  can  be  taken. 

This  is  done  by  using  the  parallel  light  from  the  top  light  path  and  blocking 
the  bottom  light  path  at  the  compensation  chamber.  The  same  optics  and  light 
sources  are  used.  Two  light  sources,  a  100  watt  continuous  mercury  arc  lamp 
with  a  d.c.  power  supply  and  a  spark  source,  can  be  mounted  on  the  interfero¬ 
meter  together  with  the  appropriate  lens  systems.  The  continuous  light  source 
is  vised  for  streak  interferograms  and  high-speed  shadowgraph  movies  while  the 
spark  light  source  is  used  for  still  interferograms  and  shadowgraph  pictures. 

The  flew  field  in  the  nozzle  is  recorded  by  shadowgraph  methods  in  par¬ 
allel  light  which  is  collimated  by  means  of  a  f/6,  610  mm  focal  length  Bausch 
&  Lomb  Aero  Tessar  aerial  camera  lens.  A  spark  light  source  with  a  0.08  cm 
diameter  source  is  located  at  the  focal  point  of  the  lens.  The  spark  source 
unit  has  a  duration  of  0.5  x  10-6  sec.  and  is  of  the  type  suggested  by 
Kovasznay*.  A  detailed  description  of  the  spark  source  and  its  power  supply 

*Kovasznay,  L.  S.  G.  1949  High  power  short  duration  spark  discharge.  Rev. 
of  Sci.  Inst.  20,  696. 


21 


A  EDC-T  R-70-7 1 


used  here  may  be  found  in  Klikoff*. 

The  spark  source  is  triggered  by  either  one  of  two  methods,  a  light 
screen  shock  detection  station  or  a  contact  on  the  diaphragm  plunger.  The 
first  method  works  in  the  following  manner.  After  rupture  of  the  diaphragm, 
the  shock  wave  travels  into  the  dump  tube  and  triggers  the  light  screen  by 
bending  a  light  beam  into  a  phototube.  The  output  of  the  light  screen  is 
then  fed  into  a  time  delay  unit  (General  Radio  Time  Delay  Generator,  Type 
1392-A)  which  in  turn  triggers  a  thyratrcn  amplifier  and  spark  source  combina¬ 
tion  after  any  desired  time  interval.  A  detailed  description  of  the  light 
screen  arrangement  and  signal  amplification  instrumentation  developed  at  our 
laboratory  may  be  found  in  Eberstein**.  The  second  method  uses  a  contact  on 
the  diaphragm  plunger  to  close  a  simple  trigger  circuit  which  sends  a  10  volt 
signal  to  the  delay  unit  which  in  turn  triggers  the  thyratron  amplifier  and 
spark  source  combination.  The  signal  is  also  sent  to  a  counter  (Beckman, 
Universal  EPUT  and  Timer  Model  7360)  which  is  connected  to  the  output  of  the 
phototube.  In  this  way,  the  plunger  contact  signal  is  referenced  to  the 
phototube  output  signal  so  that  a  consistent  timing  reference  is  used. 

The  still  interferograms  are  recorded  by  a  Leica  camera  with  a  f/4.5, 

135  mm  lens  with  a  close-up  bellows  set  at  a  magnification  ratio  of  approxi¬ 
mately  1.01.  The  image  of  the  flow  is  recorded  on  Kodak  Tri-X  Pan,  35  ntn  film 
with  an  ASA  rating  of  400.  Shadowgraphs  are  recorded  on  4  x  5  in.  type  52 
Polaroid  film,  and  in  both  photographic  systems  the  exposure  is  determined 
by  the  short  duration  (0.5  x  10“®  sec.)  spark  source.  Shadowgraph  movies  of 
the  flaw  in  the  nozzle  are  obtained  by  a  Fas  tax  16  mn  Model  WF3  high-speed 
motion  picture  camera  in  the  framing  mode  at  approximately  5000  frames  per 
second.  Streak  interferograms  were  made  by  using  this  camera  in  the  streak 
mode  (by  removal  of  rotating  prism) .  The  light  source  in  each  case  is  a  100 
watt  mercury  arc  lamp  driven  by  a  50  volt  d.c.  power  supply. 


*KLLkoff,  W.  A.  1965  Propagation  of  weak  conical  disturbances  in  relaxing 
supersonic  flows.  D.  Eng.  Thesis,  Yale  University,  New  Haven. 

**Eberstein,  I.  J.  1966  Evaluation  of  light  deflection  technique  for  detec¬ 
tion  of  transient  shock  waves.  Rev.  of  Sci.  Inst.  37,  959. 
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APPENDIX  II 
TABULATED  RESULTS 


Spark  Shadowgraph  Experiments 

Spark  shadowgraph  pictures  of  the  flew  in  the  M  =  1.67  3°  -  1/2  angle 
nozzle  without  a  model  present  were  taken  during  the  starting  process  using 
dry  air.  The  initial  pressure  in  the  high  pressure  section  (P4)  was  1  atm  at 
21.1°C.  An  initial  pressure  ratio  (P4/Pi)  of  7.6  was  used.  Pictures  were 
taken  every  .05  msec  starting  at  .29  msec  after  diaphragm  rupture  to  steady 
flew.  The  field  of  view  contained  the  nozzle  throat  and  3"  downstream.  Spark 
shadowgraph  pictures  were  also  taken  with  the  nozzle  exit  and  2"  upstream  of 
the  exit  in  the  field  of  view  by  moving  the  nozzle  4  in.  upstream.  No  shock 
waves  were  seen  during  these  experiments. 

Another  set  of  spark  shadowgraph  pictures  of  the  nozzle  exit  where  a 
.625  in.  spherical  model  was  placed  was  taken.  The  pictures  shew  the  forma¬ 
tion  of  the  steady  state  detached  shock  wave  which  forms  around  the  front  of 
the  model.  No  starting  shocks  which  pass  across  the  model  were  seen.  In 
Figure  II.  1  eight  spark  shadowgraph  pictures  of  the  flew  over  the  spherical 
model  during  the  starting  process  are  shewn.  Flew  is  from  left  to  right.  The 
first  picture  which  was  taken  at  t  =  1.85  msec  after  diaphragm  rupture  shows 
the  shear  layer  at  the  exit  of  the  nozzle.  As  time  progresses,  the  detached 
shock  is  formed  and  finally  reaches  its  steady  state  configuration  at  t  =  4.53 
msec  after  diaphragm  rupture.  The  last  picture  at  t  *  7.08  msec  is  included 
to  show  that  the  flow  is  steady  by  a  carparison  with  the  previous  picture. 

High  Speed  Shadowgraph  Movies 

High  speed  shadowgraph  movies  of  the  flew  through  the  nozzle  during  the 
starting  process  were  made  by  using  a  Fas  tax  WF-3  high  speed  16  mm  camera  at 
a  framing  rate  of  5000  frames  per  second.  Either  dry  air  or  N2  at  21.1°C  was 
used  in  the  experiments.  The  experiments  were  made  with  the  initial  pressure 
(P4)  in  the  high  pressure  tube  at  3.04  atm  and  various  pressure  ratios  (P4/P1) 
between  7.6  and  46.  Two  different  nozzle  positions  were  used  in  all  experi¬ 
ments  so  that  the  flew  between  the  nozzle  throat  and  exit  could  be  observed. 

A  similar  set  of  experiments  were  also  made  with  a  .5"  spherical  model  at  the 
nozzle  exit.  Experiments  were  performed  with  the  M  =  1.67  3°  -  1/2  angle  noz¬ 
zle  and  the  M  =  3  nozzle.  Results  of  these  experiments  shewed  no  starting 
shock  waves. 

To  improve  the  qualitative  aspects  of  the  above  experiments  the  follow¬ 
ing  modifications  were  made.  Three  1/8"  x  2"  thin  tape  strips  were  placed 
along  the  lower  M  =  1. 67  nozzle  block  and  perpendicular  to  the  flew.  The 
strips  were  placed  one  inch  apart  starting  1/2"  downstream  frem  the  throat. 
These  strips  produced  weak  disturbances  that  could  be  seen  in  our  shadowgraph 
pictures  during  the  sonic  and  supersonic  portions  of  the  starting  process. 

This  made  possible  a  very  accurate  measurement  of  starting  times.  These  ex¬ 
periments  showed  that  supersonic  flow  can  be  established  very  smoothly  and 
without  shocks  in  a  nozzle. 
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Figure  II. 1  Tine  after  diaphragm  rupture  for  above  shafkjwgraph  pictures  are; 
(a)  1.85  msec,  (b)  2.90  msec,  (c)  3.37  msec,  (d)  3.51  msec, 

(e)  3.80  msec,  (f)  4.08  msec,  (g)  4.53  msec,  (h)  7.03  msec. 

Gas:  Dry  air;  ifozzle:  M  =  1.67.  Flow  is  left  to  right. 
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Streak  Inter f erogr ams 

Streak  interferograms  of  the  starting  process  in  the  M  -  1,67  3°  -  1/2 
angle  nozzle  and  in  the  M  =  3  nozzle  were  made  by  using  the  Fas  tax  WF3  high 
speed  camera  in  the  streak  mode.  The  equivalent  film  speed  was  380  rr/sec  with 
a  time  resolution  of  approximately  5  x  10- 6  seconds.  The  experiments  were 
performed  with  dry  air  or  nitrogen  at21.1°C.  Various  initial  pressures  (P4) 
ranging  fran  1  atm  to  3.3  atm  were  used  together  with  pressure  ratios  of  46, 
15/  and  7.6.  Two  different  nozzle  positions  were  used  in  all  experiments  so 
that  interferograms  of  the  flow  between  the  nozzle  throat  and  exit  were  ob¬ 
tained.  The  quantitative  results  of  these  experiments  were  density  time  his¬ 
tories  at  various  positions  downstream  of  the  nozzle  throat  during  the  start¬ 
ing  process.  These  shew  that  the  density  undershoots  its  steady  state  value 
during  the  starting  process.  It  also  shews  that  the  undershoot  is  most  pro¬ 
nounced  at  and  near  the  throat  and  that  it  diminishes  with  increasing  down¬ 
stream  distance  until  finally  at  the  exit  it  is  not  present.  Density  profiles 
along  the  nozzle  centerline  between  the  throat  and  exit  for  a  particular  time 
during  the  starting  process  are  also  obtained  from  these  experiments.  Since 
they  do  not  contain  any  discontinuities  in  density,  no  shock  waves  were  de¬ 
tected  during  the  starting  process. 

■Static  Pressure  Measurements 

Static  pressure  measurements  of  the  supply  pressure  were  made  between 
the  exit  of  the  high  pressure  tube  and  the  transition  section.  Static  pres¬ 
sure  measurements  of  the  pressure  in  the  nozzle  were  also  made  at  three  posi¬ 
tions  in  the  nozzle;  at  the  nozzle  throat,  near  the  exit  and  just  outside  the 
exit.  The  measurements  of  the  static  supply  pressure  show  that  the  initial 
pressure  (P4)  first  undershoots  and  then  slightly  overshoots  the  steady  flow 
supply  conditions  before  it  becomes  steady.  Measurements  made  in  the  nozzle 
show  that  the  undershoot  effect  is  present  at  the  nozzle  throat  and  diminishes 
with  downstream  distance  until  finally  at  the  nozzle  exit  is  is  not  present, 
that  is,  at  the  exit  the  pressure  decreases  monotonically  to  the  steady  state 
static  pressure.  These  experiments  were  nude  using  N2  at  21.1°C  at  initial 
pressure  (P4)  ranging  from  1  atm  to  3.04  atm  with  various  pressure  ratios 
between  7.6  and  46.  In  all  cases  the  initial  pressure  and  pressure  ratios  had 
no  effect  on  the  undershoot. 

In  Figure  II. 2,  four  oscilloscope  traces  of  the  static  pressure  measure¬ 
ments  of  the  supply  tube  are  shown.  The  gas  used  in  these  experiments  was 
dry  air  at  1  atm  and  21.1°C  in  the  supply  tube.  The  pressure  in  the  lew  pres¬ 
sure  tube  was  varied  between  .079  atm  and  .289  atm  so  that  initial  pressure 
ratios  (P4/P1)  between  12.65  and  3.46  could  be  obtained.  The  vertical  scale 
is  .5  psi/div  and  the  horizontal  scale  is  2  msec/div  for  all  traces.  Results 
of  these  experiments  show  that  the  undershoot  is  not  effected  by  the  pressure 
ratio  as  long  as  supersonic  flow  is  established  in  the  nozzle.  The  steady 
state  supply  pressure  is  not  effected  by  boundary  layer  growth  in  this  tube 
as  can  be  seen  from  the  pressure  traces.  The  effect  of  the  reflected  expan¬ 
sion  wave  is  clearly  shown  in  these  pressure  traces. 

In  Figure  II .3,  three  oscilloscope  traces  of  the  static  pressure  measure¬ 
ments  of  the  supply  tube  are  shown.  The  gas  vised  was  dry  air  at  21. 1°C.  In 
the  first  trace,  the  initial  pressure  (P4)  was  3.38  atm  and  the  low  pressure 
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t  (2  msec/div) 

Figure  II. 2  Supply  pressure  vs.  time  at  upstream  position.  Initial  pressure, 
P4,  is  1  atm.  while  low  pressure,  Plf  is  varied  to  yield  different 
pressure  ratio.  Gas:  Dry  air;  Nozzle:  M  =  1.67. 
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Exp.  122-3 
Pjt/Pi  =  25.8 


Exp.  115-33 
P4/P1  =  7.6 


•H 

8. 
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Exp.  115-33a 
P4/P1  =  7.6 


t(10.0  msec/div) 


Figure  II. 3  Supply  pressure  vs.  time  at  upstream 

position.  Gas:  Dry  air;  Nozzle:  M  =  1.67 
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(Pj )  was  .131  atm  so  that  the  pressure  ratio  (P4/P1)  was  25.8.  The  undershoot 
again  appears  and  is  clearly  not  a  function  of  initial  pressure  or  pressure 
ratio  as  can  be  seen  by  ccnparing  it  with  the  second  trace  which  was  taken 
with  P4  =  1  atm  and  P4/P1  =  7.6.  The  third  trace  was  taken  with  the  gas  pres¬ 
sure  and  pressure  ratio  the  same  as  in  the  second  picture,  but  the  time  scale 
was  changed  so  that  a  long  time  history  was  recorded.  This  trace  shows  the 
starting  of  the  second  stead/  state  period.  However,  it  is  interrupted  by  the 
reflected  shock  wave  which  started  when  the  diaphragm  ruptured  and  reflected 
off  the  end  plate  of  the  lew  pressure  tube. 


28 


AE  DC-T  R-70-71 


APPENDIX  III 

DETAILS  OF  NUMERICAL  CALCULATIONS 


A  numerical  technique  was  used  to  calculate  the  gas  properties  in  the 
nozzle  during  the  starting  process  using  the  method  of  characteristics.  The 
following  assumptions  were  made  in  the  analysis;  that  the  gas  is  inviscid, 
non-conducting  and  has  constant  specific  heats  and  that  the  flew  is  adiabatic 
and  one-dimensional  with  area  change.  It  was  also  assumed  that  the  flew  of 
the  gas  in  the  high  pressure  tube  is  initiated  by  a  centered  expansion  wave 
at  the  diaphragm  position.  This  assumption  is  valid  only  if  the  characteris¬ 
tic  time  for  diaphragm  rupture  is  much  less  than  the  characteristic  time  of 
the  nozzle  flow  which  is  the  time  a  sound  wave  takes  to  travel  the  length  of 
the  nozzle.  In  this  case  the  time  for  diaphragm  rupture  is  less  than  .1  msec 
and  the  characteristic  time  of  nozzle  flow  is  1  msec.  The  calculations  were 
made  using  the  non-dimensional  characteristic  equations.  A  derivation  of  the 
characteristic  equations  and  details  of  the  numerical  technique  are  given  be¬ 
low. 

Derivation  of  Characteristic  Equations 

A  brief  derivation  of  the  characteristic  equations  follcw,  however,  for 
a  more  detailed  derivation,  see  Reference  8.  The  conservation  equations  of 
mass  and  momentum  for  one-dimensional  flow  with  area  change,  A(x) ,  are: 


ap 

3t 


+ 


U 


3p 

3X 


+  P 


3U 

3x 


+  pu 


dfcnA  _  rt 
dx  ~  U 


(1) 


3u  3u  _  1  3p 

3t  +  U  3X  p  3X 


(2) 


and  the  thermal  equation  of  state  for  a  perfect  gas  is 


p  =  pRT 


where  u  is  the  particle  velocity,  p  the  density,  A  the  area,  R  the 
gas  constant,  p  the  pressure  and  T  the  absolute  temperature.  For  the 
remaining  derivation,  subscripts  will  be  used  to  denote  partial  derivatives. 

For  isentrqpic  flew, 


fran  which 


P  =  P  (P/  sQ) 


dp 


dp  =  a2  (p)dp 


(3) 
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where  a2  is  defined  as  the  equilibrium  sound  speed. 
Using  equation  (3) ,  equation  (2)  becomes, 


u.  +  u  u  +  iiiifii 


u  v>  .  - p  =o 

t  X  p  Kx 


(3a) 


By  multiplying  equation  (1)  by  +  —  and  adding  (3a) ,  the  following  equa¬ 
tions  are  formed,  p 


ut  +  (u  +  a)ux  +  ^  [pt  +  (u  +  a)px]  =  -  ua  Un  A)x  (4) 

ut  +  (u  -  a)ux  -  J  [pt  +  (u  -  a)px]  =  ua  (Jin  A)x  (5) 

Two  characteristic  curves  are  defined  by  specifying  their  slopes.  A 
right  running  characteristic  curve,  Cjx  ,  is  defined  as  a  curve  in  the  x  -  t 
plane  with  a  slope  equal  to  dx/dt  =  u  +  a  and  a  left  running  characteristic 
curve,  Cjx  ,  is  defined  as  a  curve  in  the  x  -  t  plane  with  a  slope  equal  to 
dx/dt  =  u  -  a  .  A  change  in  any  quantity,  <t>  ,  along  Crr  or  C^r  is  writ¬ 
ten  as 


d<j>  =  [<|>t  +  (u  ±  a)  4>x]dt 


A  new  variable,  w(p)  ,  is  defined  by 


w(p)  = 


rP 


a(p) 


dp 


when  dw  =  —  dp  .  For  a  gas  with  constant  specific  heats  this  becanes. 


w<P)=^a0[Q3hri-1] 


or 


w(a)  = - T  (a  -  a  ) 

y  -  ±  o 


Substituting  this  into  equations  (4)  and  (5)  and  multiplying  each  by  dt, 
one  has, 
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[u^,  +  (u  +  aJu^dt  +  fwt  +  (u  +  a)wx]dt  =  -  ua  (£n  A)x  dt 

[ufc  +  (u  -  aju^ldt  -  Iwt  +  (u  -  a)wx]dt  =  ua  Un  A)x  dt 

By  denoting  changes  along  a  characteristic  with,  a  delta,  6  ,  the  above 
equations  are  rewritten  as, 

6  (u  +  w)  =  -  ua  (in  A)  dt  (6) 

IT  X 


6  (u  -  w)  =  ua  (Un  A)  dt  (7) 

A)  X 

where  the  subscripts  "r"  and  refer  to  the  right  and  left  characteristic 
curves.  The  equations  are  the  governing  characteristic  equations  for  isentropic 
one-dimensional  flow  with  area  change. 

Equations  (6)  and  (7)  can  be  rewritten  in  non-dimensional  form  by  using 
ao  as  a  characteristic  velocity.  Which  is  the  initial  sound  speed  of  the  gas 
in  the  high  pressure  tube,  L  as  a  characteristic  length  which  in  this  case 
will  be  the  nozzle  length,  and  L/ap  as  a  characteristic  time.  The  follow¬ 
ing  non-dimensional  variables  are  thus  defined. 


Substituting  the  above  into  equations  (6)  and  (7)  the  non-dimensional  equations 
follow: 


6  (U  +  W)  =  -  UA  (£n  A)  dx  (8) 

IT  Z 

6.  (U  -  W)  =  UA  (£n  A)  dx  (9) 

At  Z 
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The  non-dtaensional  equations  for  the  characteristic  curves  are; 


U  +  £ 


dz 

dx 


U  A 


Equations  (8)  through  (11)  and  the  appropriate  boundary  conditions  are 
all  that  is  necessary  for  a  numerical  solution. 

Numerical  Technique 

The  numerical  technique  uses  the  properties  of  a  centered  expansion  wave 
and  the  initial  known  conditions  to  calculate  the  coordinates  of  each  left 
running  characteristic  in  the  centered  expansion  wave  which  travels  into  the 
nozzle  and  also  the  gas  properties  at  each  coordinate.  The  coordinates  of  a 
point  on  a  characteristic  curve  are  determined  by  the  intersection  of  a  right 
and  left  characteristic  curve  each  starting  from  two  different  but  nearby 
points  when  the  properties  of  the  gas  are  known.  It  is  assumed  that  the  slopes 
of  the  right  and  left  characteristic  curves  do  not  change  significantly  be¬ 
tween  the  known  points  and  the  points  of  intersection.  This  assumption  is 
justified  if  the  points  are  very  near  each  other.  Since  the  characteristic 
equations  relate  the  change  in  (U  +  W)  and  (U  -  W)  along  a  left  and  right 
characteristic  curve  for  changes  in  t  ,  they  may  be  written  in  difference 
form  to  approximate  changes  in  (U  +  W)  and  (U  -  W)  for  small  At  .  This  method 
is  used  to  calculate  the  properties  of  the  gas  at  the  intersection  point  of  the 
left  and  right  characteristics  which  start  from  two  different  known  points.. 

The  two  resulting  equations  are  solved  for  U  and  W  and  hence  all  the  pro¬ 
perties  of  the  gas  at  that  point  are  kncwn.  This  process  is  repeated  until  all 
the  coordinates  of  a  particular  left  characteristic  and  the  gas  properties 
along  it  are  known.  The  next  characteristic  in  the  centered  expansion  fan  is 
then  calculated  in  a  similar  manner  by  starting  at  the  beginning  of  the  char¬ 
acteristic  and  using  the  calculated  results  along  the  previous  characteristic. 

The  details  of  the  calculation  follow.  In  Figure  III.l  a  z  -  t  diagram 
is  shown  with  the  first  left  running  characteristic  drawn.  It  is  a  straight 
line  because  the  first  left  running  characteristic  of  the  centered  expansion 
fan  is  moving  into  a  gas  which  is  at  rest  and  therefore  moves  with  a  constant 
velocity  equal  to  -  Aq  .  Along  this  characteristic  all  the  properties  of  the 
gas  are  kncwn;  they  are  the  initial  conditions  of  the  gas.  Therefore,  at  any 
point  along  this  characteristic  the  slope  of  a  right  running  characteristic 
starting  at  that  point  is  kncwn  and  is  Aq  .  In  particular,  all  the  proper¬ 
ties  at  point  2  are  known.  Since  the  flow  is  started  by  a  centered  expansion 
wave  at  z  =  0  and  i  =  0  ,  all  the  properties  of  the  gas  are  kncwn  at  point 
1  for  a  given  left  running  characteristic  starting  at  this  point  with  a  kncwn 
slope.  This  can  be  shown  by  noting  the  relation  that  exists  through  a  simple 
expansion  fan  for  the  velocity  and  sound  speed  as  related  to  the  initial  sound 
speed.  This  relation  is 
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Z-'t'DIAGRAM 


Figure  III.l 
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2 


2 

U  +  -=-=  A  = 


Y  -  1  O 


The  slope  of  the  left  running  characteristic  in  the  expansion  fan  is, 


SL  =  U  -  A 


From  the  above  equation,  U  and  A  can  be  written  in  terms  of  the  initial 
slope,  SL  ,  as  follows, 


A  = 


y  ~  1 

Y  +  1 


U  =  SL  +  A 


W  =  — — y  A 

Therefore,  by  specifying  the  initial  slope  of  the  left  running  characteristic 
which  starts  at  point  1  all  the  properties  of  the  gas  at  this  point  are  known. 

Since  the  coordinates  of  both  points  1  and  2  and  the  slope  of  the  re¬ 
spective  left  and  right  running  characteristic  are  known,  the  intersection  of 
the  two  characteristics  can  be  found  and  therefore  the  coordinates  of  that 
point.  Next,  to  find  the  properties  of  the  gas  at  point  3  the  characteristic 
equations  are  written  in  difference  form  in  the  following  way: 

along  C 

(U  +  W)  3  -  <U  +  W)  2  =  -  (UA)  2  (t3  -  t2> 

along  C£r 

(U  -  W)  3  -  (u-w)j  =  (UAh  (t3  -  tj) 

Hie  above  equations  are  then  solved  for  W3  and  U3  and  thus  all  the  pro¬ 
perties  of  the  gas  are  known  for  point  3.  The  same  procedure  is  used  to  find 
the  coordinates  of  the  next  point  on  the  No.  2  left  characteristic  and  then 
all  the  gas  properties  at  that  point.  This  is  continued  until  all  the  co¬ 
ordinates  and  gas  properties  along  this  characteristic  have  been  calculated. 
Then  the  No.  3  left  characteristic  is  calculated  using  the  same  procedure  but 
in  this  case  using  the  calculated  properties  along  No.  2  left  characteristic 
and  point  1  for  the  initial  angle  of  the  new  characteristic. 

For  the  above  calculations  the  initial  At  used  for  calculation  of  No.  2 
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left  characteristic  was  of  the  order  of  .01  and  thereafter  was  determined  by 
the  intersection  of  the  characteristics.  The  initial  slope  of  all  left  char¬ 
acteristics  starting  at  point  1  differed  frcm  the  previous  one  by  .001  degrees. 

Computer  Technique  and  Program  Flow  Chart 

The  computer  technique  consisted  of  three  basic  phases,  input,  calcula¬ 
tions  and  output.  The  input  phase  reads  in  all  pertinent  data  on  the  gas, 
such  as,  specific  heat,  initial  velocity,  and  sound  speed,  initial  increments 
in  At  and  change  in  characteristic  initial  angles,  the  nozzle  area  as  a 
function  of  z  ,  and  frequency  of  output  printouts.  The  calculations  were 
performed  for  one  characteristic  at  a  time.  After  all  data  for  a  new  char¬ 
acteristic  was  computed  it  was  printed  by  the  computer  printer  and  then  used 
as  the  kncwn  characteristic  for  the  next  calculation.  The  output  phase 
printed  the  coordinates  of  each  new  characteristic  together  with  the  oorres- 
poinding  gas  velocity,  sound  speed  and  Mach  numbers  at  that  point.  The  slope 
of  the  right  and  left  running  characteristic  was  also  printed  out  for  each  co¬ 
ordinate.  By  selection  of  the  appropriate  control  variables,  only  pre¬ 
determined  characteristics  were  printed.  A  simple  flew  chart  of  the  program 
is  shown  in  Figure  III. 2. 

A  problem  which  prevented  a  complete  solution  to  the  starting  process  was 
encountered.  This  problem  is  basically  due  to  the  divergence  of  the  left 
characteristics  near  and  at  the  nozzle  throat  when  the  Mach  number  approaches 
one.  When  this  happens  the  assumption  that  the  intersection  of  the  right  and 
left  running  characteristic  is  near  the  starting  point  of  the  two  character¬ 
istics  is  invalid.  Therefore,  the  increments  must  be  decreased  in  size  con¬ 
tinuously  to  maintain  the  validity  of  the  above  assumption.  At  this  time  only 
the  results  of  the  program  where  the  assumptions  are  valid  are  used.  Figure 
III.l  shews  the  divergence  of  the  characteristics  at  and  near  the  nozzle 
throat  as  Mach  1  is  approached. 
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Simplified  Flow  Chart 


Figure  III. 2 
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